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STAGNATION  POINT  FLOW 
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Chairman:     K.  T.  Millsaps 

Major  Department:     Aerospace  Engineering 

The  flow  at  and  in  the  neighborhood  of  a  stagnation 
point  is  obtained  from  a  known  solution  of  the  boundary 
layer  flo\\^  downstream.     Tiie  flow  is  assumed  to  impinge 
on  a  flat  plate  at  right  angles.     Tlie  solution  is  ob- 
tained using  Stokes'  approximations,  and  is  obtained  in 
terras  of  arbitrary  constants.     These  constants  are  eval- 
uated by  requiring  the  shear  stress  at  the  wall  for  the 
Stokes  flow  to  match  that  of  the  boundary  layer  flow. 

Graphs  are  obtained  by  matching  the  two  solutions 
at  1.0,  5.0,   10.0,  and  14.3  cm.   from  tlie  stagnation  point. 
A  graph  is  also  obtained  for  the  variation  of  the  shear 
stress  at  the  wall  with  distance  down  the  wall. 
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CHAPTER  I 
Introduction 


Blood  coagulation  has  been  investigated  quite  exten- 
sively by  many  people  with  volumes  of  books  written  on  the 
subject.     Although  the  process  o£  coagulation  is  a  sequence 
of  chemical  reactions,  chemistry  fails  to  give  a  complete 
answer.     The  complete  answer  appears  to  involve  hydrodynamics 
as  evidenced  by  the  recent  investigations  of  Kantrowitz.^ 
His  experiments  have  shown  that  blood  clotting  occurs  at  the 
hydrodynamic  stagnation  point. 

The  stagnation  point  on  a  body  in  a  flow  is  that  point 
at  which  the  velocities  and  the  shear  stress  are  zero.  Imme- 
diately downstream  of  the  stagnation  point  the  flows  arc 

diverging  from  the  point  (Figure  1). 
Along  the  streamline  which  inter- 
sects the  body  at  the  stagnation 


J 


'J 


^    '    ^  f  '  f      point,  the  velocity  component  par- 
Figure       1.  Flow  J^llgl  l^Q^y  .g 

The  equations  governing  the  flow  of  an  incompressible 


viscous  fluid  are  equations  derived  by  C.  L.  M.  H.  Navier 

2 

and  Sir  G.  G.  Stokes,      They  are  derived  for  a  Newtonian 
fluid  in  which  the  relationship  between  the  stress  tensor 
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and  the  rate  of  strain  tensor  is  linearized.  Although 
blood  is  not  a  Newtonian  fluid,  near  the  wall  the  blood 
particles  are  generally  absent  so  the  fluid  may  be  con- 
sidered a  plasma. 

A  solution  of  the  Navier-Stokes  equations  was  obtained 
by  Hiemenz  but  the  boundary  conditions  imposed  are  far  from 
ideal. ^    The  physical  case  would  consist  of  a  fluid  flowing 
toward  a  body,  such  as  a  flat  plate,  which  is  then  deflec- 
ted by  the  plate.     The  solution  of  the  Navier-Stokes  equa- 
tions is  subjected  to  the  boundary  condition  that  at  infin- 
ity the  velocity  component  parallel  to  the  wall  is  not  zero. 
In  other  words,  the  flow  starts  out  diverging  from  the  stag- 
nation streamline  through  no  effect  of  the  plate. 

Additional  methods  used  to  investigate  stagnation  point 
flows  involve  the  solution  of  an  asymptotic  form  of  the 
Navier-Stokes  equations.^'  ^    The  assumptions  and  the  equa- 
tions are  developed  in  Chapter  II,  and  it  is  there  stated 
that  at  the  stagnation  point  the  assumptions  of  the  theory 
are  violated, thus  rendering  the  solutions  invalid. 

Because  of  the  evidence  of  Kantrowitz  and  others,  it 
is  felt  that  a  thorough  investigation  of  stagnation  point 
flows  is  necessary.  The  problem  is  idealized  by  assuming 
a  two-dimensional  flow  to  be  impinging  at  right  angles  onto 
a  flat  plate.  The  method  to  be  used  will  be  one  put  for- 
ward by  G.  F.  Carrier  and  C.  C.  Lin.^     It  consists  of 
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obtaining  a  Stokes  flow  solution,  which  is  valid  near  the 
stagnation  point,  and  matching  this  solution  to  the  bound- 
ary layer  solution  valid  downstream  of  the  stagnation 
point.     The  assumptions  and  the  theory  are  presented  in 
Chapter  II. 

In  Chapter  III  the  solution  of  the  asymptotic  equa- 
tions valid  downstream  of  the  stagnation  point  is  pre- 
sented.    The  work  was  done  by  M.  B.  Glauert,  and  an  exact 
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solution  was  found.  The  solution  has  been  experimentally 
verified  by  P.  Bakke  with  negligible  error  present  between 
theory  and  experiment. 

Chapter  IV  is  devoted  to  the  solution  for  the  flow  in 
the  neighborhood  of  the  stagnation  point.     A  first  order 
correction  is  obtained  and  added  to  the  solution.  The 
solutions  are  obtained  in  terms  of  arbitrary  constants 
which  are  evaluated  by  matching  to  the  solution  of 
M.  B.  Glauert.     This  is  done  in  Chapter  V,  the  results 
being  contained  in  Chapter  VI, 


CHAPTER  II 
Theory  and  Method  of  Solution 

A.  Theory 

In  the  analysis  o£  an  incompressible  viscous  flow,  the 
use  of  a  dimensionless  parameter,  the  Reynolds  number,  per- 
mits the  Navier-Stokes  equations  to  be  considered  in  asymp- 
totic form,  depending  upon  the  magnitude  of  the  Reynolds 
number.  For  large  Reynolds  numbers  the  Navier-Stokes  equa- 
tions become  the  asymptotic  equations  of  L.  Prandtl's  bound- 
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ary  layer  theory,    whereas  for  small  Reynolds  numbers  the 
asymptotic  equations  describe  a  flow  called  Stokes  flow  be- 
cause they  were  originally  studied  by  Sir  George  G.  Stokes. 

The  Reynolds  number  is  the  ratio  of  the  inertial  forces 
to  the  viscous  forces.     It  may  be  found  from  the  relation 
Vd/x/  ,  where  V  is  the  velocity  of  the  fluid,  d  is  a  length 
characteristic  of  the  body,  and        is  the  kinematic  viscos- 
ity.    For  the  case  of  stagnation  point  flow,  the  body  is  a 
flat  plate  of  infinite  extent,  and  thus  has  no  fixed  length 
which  would  characterize  the  body.     As  in  flow  over  a  flat 
plate,  the  length  is  taken  as  the  distance  measured  along  the 
plate  from  the  origin  of  the  coordinate  system,  chosen  to  be 
the  stagnation  point.     Thus  in  stagnation  point  flow,  the 
Reynolds  number  is  a  "local"  one. 
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L.  Prandtl  hypothesized  that  at  large  Reynolds  numbers 
the  effects  of  viscous  forces  are  confined  to  a  narrow  re- 
gion close  to  the  body.     The  solution  of  the  flow  problem 
consists  of  a  solution  of  an  asymptotic  form  of  the  Navier- 
Stokes  equations.     This  solution  is  required  to  match  the 
inviscid  flow  past  the  body  at  a  short  distance  above  the 
body.     The  equations  are  derived  in  many  texts  and  are 
stated  here  for  future  reference."^ 

^  ="  °  ,  (2.1b) 

The  velocities  have  been  made  dimensionless  with  respect  to 

a  characteristic  velocity,  denoted  by  U,  and  the  lengths  by 

the  ratio  t)/u.     The  pressure  is  thus  dimensionless  with  re- 
2 

spect  to  pU  . 

The  solution  of  equations   (2.1a),   (2.1b),  and  (2.1c) 
for  stagnation  point  flow  by  M.  B.  Glauert  is  presented  in 
Chapter  III. 

Sir  G.  G.  Stokes  hypothesized  that  for  very  small 
Reynolds  numbers  the  inertial  forces  in  a  fluid  are  negli- 
gible when  compared  to  the  viscous  forces.     In  this  problem 
the  Reynolds  number  is  small  because  both  velocities  and  the 
physical  length  are  small.     The  method  of  obtaining  a  solu- 
tion consists  of  solving  the  Navier-Stokes  equations  in 
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which  the  inertial  terms  have  been  deleted.     Corrections  to 
this  solution  are  hopefully  obtained  by  an  iterative  tech- 
nique.    Experimental  investigations  by  D.  Tritton  show  that 

the  solution  may  be  considered  valid  in  any  region  in  which 

12 

the  Reynolds  number  is  less  than  one. 

The  Navier-Stokes  equations,  from  which  the  asymptotic 

equations  are  obtained,  are  derived  in  many  texts  and  are, 

13 

in  cartesian  coordinates 


C2.2a) 


and  the  continuity  equation  is 

The  left-hand  side,  which  contains  the  squares  of  the 
velocities,  is  the  inertial  force  on  a  fluid  element  and 
the  right-hand  side  contains  the  pressure  forces  and  the 
viscous  forces.     The     x    coordinate  is  the  distance  along 
the  plate  and  the  y  coordinate  is  the  distance  perpen- 
dicular to  the  plate  with  the  origin  at  the  stagnation 
point.     The  u  and  v  velocities  are  in  the  x  and  y  direc- 
tions respectively. 
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The  variables  can  be  made  dimensionless ,  as  in  the 
boundary  layer  equations,  by  setting 

where      is  the  kinematic  viscosity,  U  is  a  reference 
velocity,  and  ^  is  the  stream  function  defined  by 
G    ■   ^^/d  y.    V    -  -^f/^-y.. 

The  elimination  of  the  pressure  from  equations  (2.2a) 
and  (2.2b)  allows  the  Navier-Stokes  to  be  written 

A  coordinate  transformation  is  made  to  make  equation  (2.3) 
more  tractable.     The  transformation  is  to  a  complex  plane 
given  by 

z    =    X  +  ly,  2      =    X  -  ly. 
Upon  introduction  of  the  stream  function  and  noting  that 

u  +  iv    =     -2i^^  /  ^  z*,  (2.4) 
equation  (2,3)  becomes 

57^2*' "    [     ^^'^2*    at  a^^i^^j  .  (2.5) 

The  left-hand  side  now  contains  the  viscous  force  terms 
and,  using  the  notation  of  Carrier  and  Lin,  will  be  denoted 
by  L(  'f'  ),  whereas  the  right-hand  side,  which  contains  the 
inertial  force  terms,  will  be  denoted  by  L*C  *f' )  •  must 
be  noted  that  although  the  asterisk  denotes  the  complex 
conjugate  of  the  coordinate  z,  it  does  not  signify  such  in 
regards  to  the  operator  L. 


8 


B.  Method  of  Solution 

The  analytical  method  to  obtain  the  Stokes  flow  near 

14 

the  stagnation  point  is  that  put  forward  by  Stokes.  It 
consists  of  an  iterative  scheme  applied  to  equation  (2.5). 
The  first  solution  is  that  for  which  the  inertial  force 
terms  have  been  deleted,  i.e., 

L(t)=0.  (2.6a) 
Once  this  solution  is  obtained,  the  resulting  corrections 
are  obtained  by  solving  the  following  equations  in 
sequence: 

H'i>.)=  (2.6b) 

The  above  equations  are  developed  in  Chapter  IV.  Over  any 
region  in  which  the  series 

converges,  then  this     ^  will  be  taken  as  the  solution  of 
the  Stokes  flow. 

Each  of  the  equations  (2 . 6a) , , . , , (2 . 6n)  can  be  inte- 
grated to  within  an  arbitrary  function.     These  functions 
will  be  evaluated  by  requiring  that  the  solutions  meet  the 
imposed  boundary  conditions.     The  operator  L  is  the  bi- 
harmonic  operator  and, thus,  the  arbitrary  functions  must 
be  biharmonic  functions. 
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The  solutions  will  be  obtained  in  terms  of  arbitrary 
constants  which  will  be  evaluated  by  matching  to  the  bound- 
ary layer  solution.     The  method  of  matching  will  be  . 
discussed  in  Chapter  V. 


CHAPTER  III 
The  Boundary  Layer  Solution 
A.     The  Exact  Solution 

An  exact  solution  o£  the  boundary  layer  equations  was 
given  by  M.  B.  Glauert.-^^    His  work  is  summarized  in  the 
following  discussion. 

A  jet  impinging  on  a  flat  plate  is  approximated  by 
considering  a  jet  in  which  the  velocity  distribution  does 
not  vary  along  its  length.     The  jet  is  assumed  to  flow 
into  a  similar  fluid  at  rest  with  respect  to  the  plate. 
Sufficiently  far  from  the  jet  the  Reynolds  number  is  large 
enough  to  validate  a  solution  of  Prandtl's  asymptotic 
equations.     If  the  pressure  is  assumed  everywhere  con- 
stant, the  boundary  layer  equations   ( 2 . la) - C 2 . Ic)  become 

.  ^  .  (3.1a) 

,  <^ir      ^  (3.1b) 

The  boundary  conditions  to  be  satisfied  are 

CI)      )/'=u    =     v    =     Oaty    =  0 
(2)  u  — ♦  0     as    y  — ^  ^ 
A  similarity  solution  is  sought  for  which  the  velocity 
u  is  assumed  proportional  to  x*^  and  the  jet  thickness    S  is 
assumed  proportional  to  x'^ .     Substitution  of  the  assumption 
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into  equation  (3.1a)  shows  that  a  similarity  solution 
exists  provided  that 

a  +  2b    "     1.  (3.2) 
A  stream  function  is  thus  assumed, 

where 

y  =  ("■^)^X'  .  (3.3b) 
With  equations  (3.3a)  and  (3.3b)  the  momentim  equation 
(3.1a)  becomes 

where 

and  the  primes  denote  differentiation  with  respect  to    >^  . 
The  boundary  conditions  become 

/(O)     =  /  (0)     =     0  and    /(<»  )     =  0. 

Integration  of  (3.4)  between  the  limits  of     -y  and 
infinity  yields 


where  ^ 

Multiplication  of  this  equation  by    /     and  integration 
between  the  same  limits  yield  ^ 

If  no  reverse  flow  is  assumed  to  exist,  then  as    »^  decreases 

^  will  be  zero  at  any  point  ^    is  zero.     Since  the  al)0ve 
integral  will  then  be  non-zero,  the  equation  requires  that 
V  -  ;j  and  thus 

a    -     -1/2,  b     =  3/4. 
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An  exact  solution  is  found  for  equation  (3.4)  in 
terms  of  a  parameter  h.     The  solution  is 

/=  i\ 

With  the  flux  of  exterior  momentum  flux,  F,  which 
for  the  wall  jet  equals  «/*U/40,  one  can  calculate  the 
following  flow  quantities. 


(3.7c) 


B .     The  Series  Solution 

Anticipating  that  the  Stokes  flow  solution  will  be  a 
series,  equation  (3.6)  will  be  unsuitable  for  matching  on 
a  term  by  term  basis.     A  matching  of  this  type  will  be 
possible  only  for  a  series  solution  of  the  boundary  layer 
equations.     Accordingly,  assume  ^  to  be  a  power  series  of 
the  form  oo 
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Equation  (3.4)  is  readily  written  in  terms  of  the 


power  series. 


^.o  ^     ^  ^  (3.8) 

The  expansion  of  the  product  series  in  equation  (2.8) 
yields  for  the  coefficient  of 

To  have  equation  (2,8)  identically  equal  to  zero 
for  any      it  is  necessary  that  the  coefficients  of 
separately  equal  zero.     This  will  be  true  provided  that 


The  boundary  condition  of  zero  stream  function  and 
zero  velocity  at  the  wall  require  b^    =    b-j^     »  0.  The 
coefficient  b2  will  be  the  shear  stress  at  the  wall. 
Comparison  of  this  solution  with  the  exact  solution, 
equation  (2.6),  shows  that  b2     =     1/9.     With  these  three 
coefficients  all  the  rest  may  be  found. 

This  solution,  when  transformed  to  the  polar  plane, 
shows  that  the  terms  appear  as  r^©*,  r^d*^,  r^©^,...  for 
small  angles.     The  use  of  a  small  angle  approximation  was 
employed  by  Carrier  and  Lin.     They  thus  assumed  that 
although  boundary  layer  tlieory  implies  large  incrtial 
forces,  sufficiently  close  to  the  plate,  the  dominant 
forces  are  viscous. 


CHAPTER  IV 
The  Stokes  Flow 
A.     Development  of  Solution  Scheme  and  Boundary  Conditions 

As  mentioned  in  Chapter  II,  the  original  assumption  of 
Stokes  was  that  for  small  Reynolds  numbers  the  inertial 
forces  are  negligible.     This  is  not  to  say  they  are  zero; 
rather,  that  the  velocities  are  small  enough  to  neglect 
their  squares.     Assume  the  stream  function  to  be  composed 
of  additive  parts  where  the  first  term  contains  no  inertial 
forces,  the  second  contains  them  to  the  first  order,  etc., 
i.e.  , 

.  (4.1) 
Substitution  of  this  series  into  equation  (1.5)  yields 

or,  in  terms  of  the  operators  L  and  L  , 

Since  the  function    ^  is  assumed  to  contain  no  iner- 
tial effects,   it  will  be  a  solution  of  the  equation 

/Cf^J=0.  ■  (4.2) 

This  is  the  biharmonic  equation  making  a  bihar- 

monic  function.     The  remaining  equation  to  be  solved  is 
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H  ....       =  azt). 

The  function    /l   may  be  eliminated  from  L(  ^  )  by 
setting 

equal  to  all  the  terms  of  L  ( 
which  are  functions  of    ^  only.     It  is  obvious  that  the 
equation  in  this  case  is 

^^may  similarly  be  eliminated  by  setting 
equal  to  the  remaining  terms  of  L*(  ^)  which  are  functions 
of  and  i.e., 

Addition  and  subtraction  of  the  terms 

permits  writing  equation  (4. 4a)  in  terms  of  the  operators; 

(4.4b) 

Continuing  in  like  manner  the  n-th  part  of  the  stream  func- 
tion will  be  found  from  the  equation 
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This  method  will  yield  the  solution 

over  any  region  in  which  the  series  is  convergent  and 
satisfies  the  original  assumptions  of  the  theory. 

A  solution  for  each  of  the  parts  of  the  stream  func- 
tion can  be  found  to  within  an  arbitrary  biharmonic  func- 
tion.    Denoting  the  particular  integral  by         and  the 
arbitrary  function  by  ,  the  solution  of  any  of  the 

equations  may  be  written 

The  coordinate  systems  to  be  used  to  obtain  the 
solutions  of  the  equations  is  the  z,2*  plane,  while  the 
boundary  conditions  are  more  readily  applied  in  the  polar 
plane,  with  the  solution  plotted  in  the  cartesian  plane 
(Figure  2) . 


Figure  2,     Coordinate  System 
With  the  stagnation  point  located  at  Co»0)i  ^^e 

boundary  conditions  are 

(1)  id-  =     0  at  S»  O.TT 

(2)  u  +  iv    =    -Ix^^^l  ^  z*  =    0  at  e=c,rr 

(3)  u    =     ReC-2i  ^^^'1  3  z*)  =     0  at  tf* 
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B.     Solution  With  No  Inertial  Forces 

The  first  approximation  to  the  flow  is  given  by  the 
solution  of  equation  (4.2).     E.  Goursat^^  gave  a  solution 
of  the  equation  in  terms  of  two  arbitrary  functions  ^^(z) 
and    y  (z)  and  their  conjugates     ^*(z)  and  X*^^^' 
was  found  to  be 

which  can  be  written  in  a  more  convenient  form,  viz., 

The  series  solution  of  the  boundary  layer  equation 
resulted  in  the  terms  being  multiples  of  the  fourth  root 
of  the  independent  variables  r  and  0  .     Accordingly,  the 
function  ^iz)   is  assumed  to  vary  in  like  manner.,  A 
general  expression  can  be  written  for  ^(.z)   in  terms  of 
the  infinite  series 

Upon  substitution  of  this  series  into  equation  (4.5),  it 
is  found  that  the  first  boundary  condition  will  be  satis- 
fied for  ft 

The  general  solution  may  thus  be  written 
'0  ^-o 
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The  velocities  are  found  using  equation  (2.4)  which 
states  that  u  +  iv    =    -2i  <^  ^/  ^  z* .     Upon  differen- 
tiation and  transformation  to  the  polar  plane,  they  may 
be  found  by  the  relation 

The  boundary  condition  requiring  zero  velocities  at  ff=0 
is  identically  satisfied,  but  at  S^TT 

At    ^M^thQ  boundary  condition  can  be  satisfied  if  Cn+l)/4 
=    m,  m    =     0,1,2,....     The  third  boundary  condition  pro- 
duces a  similar  situation  in  that  at 

As  before,  this  boundary  condition  will  be  satisfied  pro- 
vided that  (m-2)/2     =  0,1,2  

A  solution  to  the  Navier-Stokes  eciuations,  in  which 
the  inertial  forces  have  been  neglected  and  which  satis- 
fies the  boundary  conditions  of  stagnation  point  flow,  is 

Jso  (4.7a)  ■ 


with 


•  ( Cos.  2£0  -  ^/AJ  /J a)  -  s/^  ^^^J^* 


(4.7b) 


and  ^  (4.7c) 
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C.    The  First  Iteration 

A  solution  is  sought  for  equation  (4.3)  subject  to  the 
mathematical  and  boundary  conditions  stated  in  Part  A  of 
this  chapter.     The  particular  integral,   j^,^  ,  will  be  ob- 
tained and  the  biharmonic  functions,  /(y  ,  will  be  eval- 
uated so  as  to  make  the  boundary  values  correct. 

Substitution  of  equation  (4.3a)  into  equation  (4.3) 
results  in 


a 


The  integration  is  readily  performed  and  the  particular 
integral,  ,  when  transformed  to  the  polar  plane,  becomes 


(4.8) 
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This  particular  integral  satisfies  boundary  condi- 
tions (1)  and  (3);  that  is,  zero  mass  flow  through  the 
wall  and  u    -    0  at  0=^/2,  but  at  the  wall 

with 

fu  +  iv)  =     -  (u  +  iv) 

Thus,  the  arbitrary  function  must  be  evaluated  to  make  the 
velocities  zero  at  the  wall. 

The  general  form  for  a  biharmonic  function  was  given 
in  equation  (4.5).    Three  possibilities  exist  for  satis- 
fying the  boundary  condition  (2) : 

CD  =    0.  ?^  0, 

C2)  =     0,     X(z)     f  0, 

(3)       <^(z).    Xiz-)     f  0. 

Case  CI):     'X(^)  ^  o 

Assume  the  biharmonic  function,    ^    ,  to  be  given  by 

It  is  readily  found  that 

with 

(u  +  iv)  =  -  (u  +  iv)  , 


Thus  ,  when  added  to  ,  will  result  in  zero  veloc- 

ities at  the  wall  provided 

Case  (2):    ^(z)     =  0 

The  biharmonic  function,  denoted  by  ^  ,  may  be 
written 

The  velocities  at  the  wall  are  found  to  be 
with 

(u  +  iv)  =       - (u  +  iv) 

This  function  will  adjust  the  value  of  ^  to  satisfy  bound- 
ary condition  (2)   at  the  wall  if 

Case  C3):     ^(z)  ,     'J^  (z)     ?^  0 

This  biharmonic  function,   j^j   ,   is  essentially  a 
linear  combination  of  cases  (1)  and  (2),  with  the  differ- 
ence being  in  the  constant. 


The  velocities  at  the  wall  arc  given  by 
with 

fu  +  iv)  =        -(u  +  IV) 

As  before,  the  solution  will  sum  to  zero  provided  that 

Each  of  the  biharmonic  functions,         ,  j  =  1,2,3, 
satisfy  the  boundary  conditions  (1)  and  (3)  and  thus 
produce  no  change  in  their  values. 

A  solution  for  stagnation  point  flow,  for  which  a 
first  order  correction  of  the  inertial  effects  has  been 
made,  thus  becomes 


'  lA-Jm"'^'""''"-'-'" 
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where  ot>        o*  ,  . 


or 

with  Aj,^  Bj^j^    and  Cj,^^  as  previously  defined. 

The  solution  of  equation  (4.4b)   results  in  terms  begin- 
ning with  the  ninth  power  of  r  with  the  coefficient  being 

const,   (times)  aj^aja^  (10  ). 
Since  the  series  is  presumed  convergent,  it  is  felt  that 
these  terms  will  be  negligible  compared  to  those  already 
calculated  and,  thus,  the  iterative  scheme  is  terminated 
at  this  point. 


CHAPTER  V 
The  Matching  Procedure 

The  arbitrary  constants  obtained  in  Chapter  IV  must 

be  evaluated  subject  to  some  boundary  condition.     The  idea 

to  be  employed  here  is  that  used  by  Carrier  and  Lin  and 
6  17 

Imai.  '     Rather  than  require  the  solution  to  have  some 
velocity  distribution  at  infinity,  the  flow  will  be  re- 
quired to  "carry  over"  into  the  boundary  layer  solution. 
This  method  of  direct  comparison  of  two  asymptotic 
solutions  is  subject  to  a  large  number  of  questions  for 
which  there  are  very  few  answers.     In  the  final  analysis, 
experimentation  must  provide  the  answer  as  to  where  this 
"carry  over"  can  occur,  if  such  a  region  does  indeed  exist. 
Some  of  the  terms  of  the  Carrier  and  Lin  solution  matched 
Blasius'  boundary  layer  solution  when  the  small  angle 
approximation  was  made.     In  this  problem,  however,  tlie 
term  containing  the  sixth  power  of  r  has  the  factor 

compared  to     $  . 

1 8 

Van  Dyke      states  several  requisites  for  matching  dif- 
ferent parameter  expansions.     Examples  are,  that  one  of  tlie 
solutions  be  convergent  over  the  whole  region,  or  that  some 
similarity  rule  be  followed.     This  solution  does  not  sat- 
isfy either  of  them,  as  the  Stokes  flow  solution  is  not 
complete,  and  a  similarity  rule  is  not  applicable. 

24 


25 


It  may  be  inferred  from  the  discussion  by  Van  Dyke 
that  unless  the  above  are  satisfied,  the  Stokes  flow  solu- 
tion should  be  matched  to  some  intermediate  solution.  An 
intermediate  solution,  between  Stokes  flow  and  boundary 
layer  flow,  would  be  a  solution  either  of  the  equations 
of  Oseen^^  or  the  full  Navier-Stokes  equations.  Oseen 
flow  is  essentially  a  perturbation  of  Stokes  flow;  and, 
for  this  problem,  was  found  to  be  of  no  use.     It  should  be 
remembered  that  there  is  a  solution  of  the  Navier-Stokes 
equations  for  this  type  of  problem,  but  the  boundary  con- 
ditions at  infinity  are  not  consistent  with  the  particular 
physical  situation. 

The  method  to  be  used  here  is  similar  to  that  used  by 
Imai.^^    lie  matched  the  skin  friction  at  a  Reynolds  number 
of  one  for  the  flow  over  a  flat  plate.     For  the  present 
problem  a  fixed  location  cannot  be  defined  at  wliicli  tlic 
Reynolds  number  is  one. 

In  the  experimental  investigation  of  M,   B.  Clauert's 
solution,  it  was  found  that  at  14.3  cm.   from  the  stagna- 
tion point,  the  difference  between  theory  and  experiment 
was  negligible . There  is  no  mention  of  measurements 
closer  than  14.3  cm.,  so  the  error  introduced  by  matching 
at  smaller  distances  from  the  origin  cannot  be  determined. 
Other  points  will  be  picked  at  which  to  match  in  order 
that  the  effect  on  location  can  be  ascertained. 
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For  each  constant  which  is  retained  it  is  necessary 
to  differentiate  the  stream  function.     Since  the  solution 
of  equation  (4.4b)  would  begin  with  terms  beginning  at  r^, 
the  series  retained  must  terminate  with  r^ ,  the  r^  term 
being  incomplete.     This  terminated  series  will  contain 
three  arbitrary  constants,  a^,  i    =    0,1,2.     Since  the 
zeroth  and  first  derivatives  are  identically  zero  at  the 
plate,  the  second  through  the  fourth  derivatives  are 
necessary. 

The  partial  operators,  which  when  evaluated  at  the 
wall,  become 


The  truncated  series  is 


^  =  /•  Y G-coisej-f-  2a,  Cos  sfi  -  Cos  s-ff) 
^  Za^r'^/'ecsJTs  -  cos  re)  ^ /^//*  ^s/aj^^  ^{^  s/4j44 


where 
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and  . 

Differentiation  of  the  series  solution  for  the  bound- 
ary layer  flow  shows  that 


and 


where  \>2  is  found  to  be  1/9. 

The  values  of  the  constants  are  contained  in  the  next 
chapter. 


CHAPTER  VI 


Results  and  Conclusions 


The  locations  of  the  matching  points  are  determined 

by  considering  the  properties  of  blood  and  the  experimental 

evidence  of  liakke.     Blood  has  been  found  to  have  a  kine- 

7  2 

matic  viscosity  of  3.8(10'^)  cm.     /sec.  with  its  average 
velocity  in  blood  vessels  to  be  in  the  range  of  30.0  to 
40.0  cm. /sec. Choosing  for  convenience  the  reference 
velocity,  U,  to  be  38.0  cm. /sec,  the  non-dimensional- 
izing  factor  for  the  lengths  becomes  1(10"-^)  cm.  Mith 
the  experimental  verification  of  the  boundary  layer  solu- 
tion at  14.3  cm.,  the  first  matching  point  chosen  is 
14,300.     Other  points  chosen  correspond  to  1.0,  5,  0,  and 
10.0  cm.     The  solutions  are  plotted  in  Figures  3-6.  Fig- 
ure 7  shows  tlie  solution  for  matching  at  14.3  cm.  on  a 
diminished  scale.     The  values  for  the  coefficients  are 
contained  in  Table  1. 

The  solutions  are  plotted  only  for  the  first  case, 


ently  about  93.7  percent  of  J^^  and  94.7  percent  of  //S 
for  tlie  matching  at  14.3  cm.     For  the  other  matching 


0.     Examination  of  the  data  shows  that 


in  the  region  about  the  stagnation  point. 


points 


is  never  less  than  92.6  percent 
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percent  of  f^j  in  the  region  of  interest.     Differences  in 
the  graphs  would  not  be  detectable  except  on  a  greatly 
expanded  scale .  ' 

Lines  of  constant  Reynolds  numbers  are  also  plotted 
to  determine  the  region  of  validity  for  the  Stokes  flow 
solution.     The  Reynolds  number  can  be  viewed  as  the  ratio 
of  the  inertial  forces  to  the  viscous  forces;  thus,  be- 
cause of  the  assumption  that  inertial  forces  are  negli- 
gible in  Stokes  flow,  the  region  of  validity  is  that  re- 
gion in  which  the  Reynolds  number  is  less  than  one.  It 
appears  at  first  that  for  very  small  distances  from  the 
origin,  the  Reynolds  number  would  be  small  regardless  of 
the  velocities.     The  assumption  of  negligible  squares  of 
the  velocities  would  be  violated  even  though  the  Reynolds 
number  is  small.     A  choice  of  the  vertical  distance  rather 
than  the  horizontal  distance  as  the  characteristic  length 
might  be  an  alternative,  but,  for  the  solutions  obtained, 
the  regions  graphed  are  all  at  a  Reynolds  number  less  than 
one.     The  choice  of  the  horizontal  distance  is  more 
restrictive  in  defining  the  region. 

If  the  limit  of  the  validity  of  Stokes  flow  is  chosen 
at  a  Reynolds  number  of  .01,  the  region  for  matching  at 
1.0  cm.   is  essentially  a  semicircle  of  radius  012  cm., 
while  for  the  matching  at  14.3  cm.,  the  region  is  essen- 
tially a  semicircle  of  radius  1.4-1.6  cm.     The  reason  is 
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quite  apparent  when  the  boundary  layer  solution  is  con- 
sidered.    The  solutions  were  matched  by  requiring  the 
shear  stress  and  its  derivatives  to  be  equal  at  some 
point.     The  shear  stress  in  the  boundary  layer  flow  is 
proportional  to  r'^/^.     Analytically,  the  shear  stress 
is  directly  proportional  to  the  velocity  gradient,  so 
that  as  the  shear  stress  decreases,  near  the  wall  there 
is  less  mass  flow. 

Beyond  a  Reynolds  number  of  about  unity,  the  flow 
begins  to  depart  from  that  which  had  been  expected;  but 
this  region  is  certainly  well  beyond  that  at  which  a 
Stokes  flow  solution  could  be  considered  valid. 

The  assumption  of  Chapter  IV,  that  the  second  iter- 
"ation  is  negligible,  is  borne  out  by  examining  Table  1. 
In  Chapter  IV  it  was  stated  that  the  first  term  of  this 
second  iteration  contained  a  triple  product  of  the  con- 
stants.  Even  if  the  product  is  aQ ,  the  term  will  be  of 
order  10~^^,  which  can  be  considered  negligible  compared 
to  aQa-j^. 

Figure  8  shows  the  variation  of  the  shear  stress  at 
the  wall  with  the  distance  down  the  wall  for  the  matching 

at  14.3  cm.     The  shear  stress  has  been  made  dimensionless 

2 

with  respect  to        U  .     The  variation  is  shown  for  the 
two  cases  of     (f  ii)     =     0  and   yl^Cz)     =     0-     The  third  case 
falls  between  these  values.     The  variation  is  essentially 
linear  with  both  cases  having  values  of  8.88(10"^)  at  the 
match  point. 
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The  stagnation  point  flow  can  thus  be  considered  as 
Stokes  flow  up  to  the  matching  location  and  then  boundary- 
layer  flow.     Figure  6  gives  the  flow  which  will  produce 
the  experimentally  verified  boundary  layer  flow. 
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